EXTREMAL DOMAINS FOR THE FIRST EIGENVALUE IN A 
GENERAL COMPACT RIEMANNIAN MANIFOLD 



ERWANN DELAY AND PIERALBERTO SICBALDI 

Abstract. We prove the existence of extremal domains with small prescribed volume for 
the first eigenvalue of the Laplace-Beltrami operator in any compact Riemannian manifold. 
This result generalizes a results of F. Pacard and the second author where the existence 
of a nondegenerate critical point of the scalar curvature of the Riemannian manifold was 
required. 
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1. Introduction and statement of the result 

Let (M, g) be an n-dimensional Riemannian manifold, Vt a connected and open domain in 
M with smooth boundary, and An > the first eigenvalue of the Laplace-Beltrami operator 
— in Vt with Dirichlet boundary condition. The domain Vt is said to be extremal (for 
the first eigenvalue of the Laplace-Beltrami operator under Dirichlet boundary condition) 
if it is a critical point for the functional Vt \ — )■ in the class of domains with the same 
volume. 

An extremal domain is characterized by the fact that the first eigenfunction of the 
Laplace-Beltrami operator with Dirichlet boundary condition has constant Neumann data 
at the boundary. This result has been proved in the Euclidean space by P.R. Garabedian 
and M. Schiffer in 1953 [1], and in a general Riemannian manifold by A. El Soufi and S. 
Bias in 2007 [2]. Extremal domains are then domains where the elliptic overdetermined 
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AgU + A M = in f2 

u > in Q 

u = on d^l 

gCVu, v) = constant on dQ 

can be solved for some positive constant A, where u denotes the outward unit normal vector 
about dQ for the metric g. 

In the only extremal domains are balls. This is a consequence of a very well known 
result of J. Serrin: if there exists a solution u to the overdetermined elliptic problem 

= in n 

> in n 

= on dn 

= constant on dQ , 

for a given bounded domain Q C M" and a given Lipschitz function /, where u denotes 
the outward unit normal vector about dQ and (■,■) the scalar product in M", then Q 
must be a ball, [19]. In the Euclidean space, round balls are in fact not only extremal 
domains, but also minimizers for the first eigenvalue of the Laplacian with Dirichlet 
boundary condition in the class of domains with the same volume. This follows from the 
Faber-Krahn inequality, 

(3) \n > ^B"{n) 

where B'^{Q) is a ball of M" with the same volume as Q, because equality holds in ([3]) if 
and only if Q = see [3] and [9]. 

Nevertheless, very few results are known about extremal domains in a Riemannian man- 
ifold. The result of J. Serrin, based on the moving plane argument introduced by A. D. 
Alexandrof in [1], uses strongly the symmetry of the Euclidean space, and naturally it 
fails in other geometries. The classification of extremal domains is then achieved in the 
Euclidean space, but it is completely open in a general Riemannian manifold. 

For small volumes, a method to build new examples of extremal domains in some Rie- 
mannian manifolds has been developed in [12] by F. Pacard and P. Sicbaldi. They proved 
that when the Riemannian manifold has a nondegenerate critical point of the scalar cur- 
vature, then it is possible to build extremal domains of any given volume small enough, 
and such domains are close to geodesic balls centered at the nondegenerate critical point 
of the scalar curvature. The method fails if the Riemannian method does not have a 
nondegenerate critical point of the scalar curvature. 

In this paper we improve the result of F. Pacard and P. Sicbaldi by eliminating the 
hypothesis of the existence of a nondegenerate critical point for the scalar curvature. In 
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An + f{u) 
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particular, we are able to build extremal domains of small volume in every compact Rie- 
mannian manifold. 

For e > 0, we denote by B^{p) C M the geodesic ball of center p E M and radius e. We 
denote by C M" the Euclidean ball of radius e centered at the origin. The main result 
of the paper is the following: 

Theorem 1.1. Let M be a Riemannian manifold of dimension n > 2. There exist eo > 
and a smooth function 

$ : M X (0, eo) — > M 

such that: 

(1) For all e G (0,eo), if p is a critical point of the function $(-,e) then there exists 
an extremal domain G M whose volume is equal to the Euclidean volume of B^. 
Moreover, there exists c > and, for all e G (0, eo), the boundary ofQ^ is a normal 
graph over dB^{p) for some function v{p, e) with 

\Hp,e)\\c2.c(^dBl{p)) < ce^ 

(2) There exists a function r defined on M that can be written as 

Y = Ki \\Riemf + K2 || fficf + K:iR^ + AgR 

where Riem, Ric, R denote respectively the Riemann curvature tensor, the Ricci 
curvature tensor and the scalar curvature of {M,g), and Ki, K2, and K4 are 
constants depending only on n, such that for all k > 

||$(p,e) - Rp- e^rpllcfc(Af) < Cfce'^ 

for some constant Ck > which does not depend on e G (0, eo) (the subscript p 
means that we evaluate the function at p). 

(3) The following expansion holds: 

. X -2 n{n + 2) + 2X1 

bn[n + /) 

where Ai is the first Dirichlet eigenvalue of the unit Euclidean ball. 

The explicit computation of the constants Ki is given in section [71 We remark that 
if M is compact, then there exists always a critical point of $(-,e), and then we have 
small extremal domains obtained as perturbation of small geodesic balls in every compact 
Riemannian manifold without boundary. 

It is clear that this theorem generalizes the result of [12] because the construction of 
extremal domains does not require the existence of a nondegenerate critical point of the 
scalar curvature. In fact, if the scalar curvature function R has a nondegenerate critical 
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point po, then for all e small enough there exists a critical point p = p{e) of $(-,e) such 
that 

dist(p,po) < ce^. 

and then the geodesic ball B^{p) can be perturbed in order to obtain an extremal domain. 
We recover in this case the result of [12] , but with a better estimation of the distance of p 
to po (in [12] the distance between p and po is bounded by ce). In particular, we have the 
p- independent expansion 

The result of [12] can not be applied to some natural metrics such that a Einstein one, i.e 
when Ric = kg for some constant k, or simply a constant scalar curvature one. In the 
case where i? is a constant function, one gets the existence of extremal domains close to 
any nondegenerate critical point of the function r. In the particular case where the metric 
g is Einstein we obtain extremal domains close to any nondegenerate critical point of the 
function (we will see that Ki ^ 0) 

p — )■ II ffiemp 11^ . 

In order to put the result in perspective let us digress slightly. The solutions of the 
isoperimetric problem 

:= min VoL. dQ 

ncM : Voig n=K ^''^ 

are (where they are smooth enough) constant mean curvature hypersurfaces (here gi^ de- 
notes the induced metric on the boundary of Vt). In fact, constant mean curvature are the 
critical points of the area functional 

fi^VoL. dVL 

under a volume constraint YoXgVL = k. Now, it is well known (see [3], [n] and [TU]) that 
the determination of the isoperimetric profile is related to the Faber-Krahn profile, 
where one looks for the least value of the first eigenvalue of the Laplace-Beltrami operator 
amongst domains with prescribed volume 

FKi^ := min Xn 

QCM : Volg n=K 

A smooth solution to this minimizing problem is an extremal domain, and in fact extremal 
domains are the critical points of the functional 

under a volume constraint Vol^ Q = n. 

The result of F. Pacard and P. Sicbaldi [12] had been inspired by some parallel results 
on the existence of constant mean curvature hypersurfaces in a Riemannian manifold M. 
In fact, R. Ye built in [22] constant mean curvature topological spheres which are close 
to geodesic spheres of small radius centered at a nondegenerate critical point of the scalar 
curvature, and the result of F. Pacard and P. Sicbaldi can be considered the parallel 
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of the result of R. Ye in the context of extremal domains. The method used in [12] is 
based on the study of the operator that to a domain associates the Neumann value of its 
first eigenfunction, which is a nonlocal first order elliptic operator. This represents a big 
difference with respect to the result of R. Ye, where the operator to study was a local 
second order elliptic operator. 

In a recent paper, [13], F. Pacard and X. Xu generalize the result of R. Ye by eliminating 
the hypothesis of the existence of a nondegenerate critical point of the scalar curvature 
function. For every e small enough, they are able to build a small topological sphere of 
constant mean curvature equal to by perturbing a small geodesic ball centered at 
a critical point of a certain function defined on M which is close to the scalar curvature 
function. For this, they use the variational characterization of constant Hq mean curvature 
hypersurfaces as critical points of the functional 

in the class of topological sphere, where -D5 is the domain enclosed by S", see [T3] . 

Our construction is based on some ideas of [IS]. For this, we use the variational char- 
acterization of extremal domains. The main difference and difficulties with respect to the 
result of F. Pacard and X. Xu arise in the fact that there does not exist an explicit formu- 
lation to compute the first eigenvalue of a domain while there exists an explicit formulation 
to compute the volume of a surface. 

Our result shows once more the similarity between constant mean curvature hypersur- 
faces and extremal domains. The deep link between such two objects has been underlined 
also in [16] and [T7] . 

It is important to remark that P. Sicbaldi was able to build extremal domains of big 
volume in some compact Riemannian manifold without boundary by perturbing the com- 
plement of a small geodesic ball centered at a nondegenerate critical point of the scalar 
curvature function, see [2U]. As in the case of small volume domains, the existence of a 
nondegenerate critical point of the scalar curvature function is mandatory (and such result 
requires also that the dimension of the manifold is at least 4). It would be interesting to 
adapt our result in order to build extremal domains of big volume in any compact Rie- 
mannian manifold without boundary by perturbing the complement of small geodesic balls 
of radius e centered at a critical point of the function $(-,e) or some other similar func- 
tion. This result would allow for example to obtain extremal domains Vt^ that are given 
by the complement of a small topological ball in a flat 2-dimensional torus, and by the 
characterization of extremal domains this would lead to a nontrivial solution of ([2]), with 
fit) = At, in the universal covering VL^ of fi^, which is a nontrivial unbounded domain of 
M^. Up to our knowledge the existence of this unbounded domain is not known. Remark 
that VL^ is a double periodic domain, made by the complement of a infinitely countable 
union of topological balls. The existence of Vt^ would establish once more the strong link 
between extremal domains and constant mean curvature surfaces, via the double periodic 
constant mean curvature surfaces (see [6], [15] and [H]). 
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2. Notations and preliminaries 

Let Qq be a smooth bounded domain in M. We say that {Qt}te{-to,to) is a deformation of 
Oo if there exists a vector field H such that Qt = C,{t, Qq) where $,{t, •) is the flow associated 
to S, namely 

^(t,p) = S(e(t,p)) and mp)=P- 

In this case we say that S is the vector field that generates the deformation. The de- 
formation is said to be volume preserving if the volume of Qt does not depend on t. If 
{Qt}t€{-to,to) is ^ deformation of Qq, and Xn^ and Ut are respectively the first eigenvalue 
and the first eigenfunction (normalized to be positive and have L'^{Qt) norm equal to 1) of 
— Ag on Qt with Dirichlet boundary condition, both applications 1 1 — > Xn^ and 1 1 — y ut 
inherit the regularity of the deformation of Qq. These facts are standard and follow at once 
from the implicit function theorem together with the fact that the least eigenvalue of the 
Laplace-Beltrami operator with Dirichlet boundary condition is simple. 

A domain is an extremal domain (for the first eigenvalue of — Ag with Dirichlet 
boundary condition) if for any volume preserving deformation {Qt}t£{-to,to) of ^o? we have 



dAo, 
dt 



Assume that {^lt}t is a perturbation of a domain generated by the vector field H. 
The outward unit normal vector field to dQt is denoted by I't- We have the following result, 
whose proof can be found in [2] or in [T2j: 

Proposition 2.1. (Garahedian - Schiffer, El Soufi - Bias). The derivative of the first 
eigenvalue with respect to the deformation of the domain is given by 

t=0 J dQo 



dXn^ 
dt 
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This result allows to characterize extremal domains as the domains where there exists a 
positive solution to the overdetermined elliptic problem 



(4) 



AgU + Xu = Q in O 



u = on dQ 



gCVu, v) = constant on dQ 



for a positive constant A, where u is the outward unit normal vector about dQ. The proof 
of this fact follows directly from Proposition I2.H but can be found also in [12] . 

Given a point p E M we denote hj Ei, . . . , an orthonormal basis of the tangent plane 
TpM. Geodesic normal coordinates x := {x^ , . . . , x") G at p are defined by 

X{x):= Exp^ ^XJ^^ 

where Exp^ is the exponential map at p for the metric g. 

It will be convenient to identify M" with TpM and S^~^ with the unit sphere in TpM. If 
X := (x^, . . . , x") G M", we set 

n 

(5) Q{x) ■=^x' EiETpM . 

i=l 

It corresponds to the vector of TpM whose geodesic normal coodinates are x. Given a 
continuous function / : S*""^ i — ?■ (0, +oo) whose L°°-norm is small (say less than the cut 
locus of p) we define 

:= |Exp^(e(x)) : x G M" < |x| < / 

For notational convenience, given a continuous function / : 5""^ — )■ (0, oo), we set 

Bf := {x G M" : < |x| < /(x/|x|)} . 

When we do not indicate the metric as a superscript, we understand that we are using 
the Euclidean one. Similarly, we denote by Vol^ the volume in the metric g, by dvol^ the 
volume element in the metric g to integrate over a domain, by dvol^j^^ the volume element 
in the induced metric gin to integrate over the boundary of a domain. When we do not 
indicate anything we understand that we are considering the Euclidean volume, or the 
Euclidean measure, or the measure induced by the Euclidean one on boundaries. 

Our aim is to show that, for all e > small enough, we can find a point p G M and a 
function v : S*""^ — > M such that 

Vol, (p) = Vol B, = Vol Bi = e-^ 
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(where Un is the Euchdean volume of the unit sphere S"- ^) and the overdetermined problem 

' A,0 + A0 = in 
(6) I = on 

^ g(y(t>,u) = constant on dBl^^^^.^{p) 

has a non trivial positive solution for some positive constant A, where v is the unit normal 
vector field about dB^^^^^^ij)). 

Clearly, this problem does not make sense when e = 0. In order to bypass this problem, 
we observe that, considering the dilated metric g := (7, the above problem is equivalent 
to finding a point p E M and a function v : 5""^ — > R such that 

Vols Bl^{p)= Vol B, 

and for which the overdetermined problem 

' Ag0 + A0 = in Bf^,{p) 
4> = on dBl^M 
giy^cp^u) = constant on 9i?f_^^(p) 



(7) 



has a non trivial positive solution for some positive constant A, where u is the unit normal 
vector field about dBf_^^{p). Taking in account that the functions (j) and (p have L^-norm 
equal to 1, we have that the relation between the solutions of the two problems is simply 
given by 

A _ ,-n/2 7 



and 



3. Some expansions in normal geodesic coordinates 

We precise that through this paper we consider the following definition of the Riemann 
curvature tensor: 

Riem{X, Y)Z = Vx^yZ - VyVxZ - V[x,y]Z 

where V denotes the Levi-Civita connection on the manifold M. 

Geodesic normal coordinates are very useful because there exists a well known formula 
for the expansion of the coefficients of a metric near the center of such coordinates, see 



EXTREMAL DOMAINS FOR THE FIRST EIGENVALUE 

21], [11] or [18]. At the point of coordinate x, the following expansion holdS 



Qij ^ij 2 ^ ^ Rikjl XX g ^ ^ Rikjl,m X XX 

k,£ k,£,m 

1 - 2 - 

— ^ ^ Rikjl,mc7 X X X X -\- -— ^ ^ Rikjt Rirnja X XX X -\- (D(\x\ ') 



20 ^ 45 

k,£,m,a k,£,m,a 



where 



Rikji = g{Rie'mp{Ek, Ei) Ej, Ei) 
Rikje,m = g{{VE^Riem)p{Ek,Ei) Ej,Ei) 

Rikji g{{VE^^E^Riem)p{Ek, Ei) Ej, Ee) , 

and the subscript p means that we evaluate the quantity at p. In ([8]) the Einstein notation 
is used (i.e., we do a summation on every index appearing up and down). Such notation 
will be always used through this paper. 

This expansion allows to obtain other expansions, as those of the volume of a geodesic 
ball, or the first eigenvalue and the first eigenfunction on a geodesic ball. In order to recall 
such expansions, let us introduce some notations. Let us denote by Ai the first eigenvalue 
of the Laplacian in the unit ball Bi with Dirichlet boundary condition. We denote by 
(j)i the associated eigenfunction 



(9) 



A01 + Ai 01 = in Bi 
01 = on dBi 



normalized to be positive and have L^(i?i) norm equal to 1. It is clear that 0i is a radial 
function 0i(x) = 0i(|x|). We denote r = 

We recall now some expansions we will need later, whose proofs can be deduced from 
([8]). We remind to [13] and [8] for the proofs. For the volume of a geodesic ball of radius e 
we have: 

(10) e-'^ Vol,5f (p) = ^ + WQe^ + We^ + O(e^), 



where 



(11) 



^0 = -7T-7—l^Rp 



Qn {n + 2] 



-"^We choose the convention of [H] , some sign in the development are different from those in [T3] or 
because of a different choice of the definition of Rtjki 
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For the first eigenvalue of tlie Laplace-Beltrami operator witli Diriclilet boundary condi- 
tion on a geodesic ball of radius e we have: 

(12) A5.(,) = Ai + Ao + A + 0(e^) 
where 

Ao - 

(13) 2 

A = (3||i?zempf + ^||i?zCp||2 + ^!^i?2 + l(A,i?)J 

n[n + 2) \ 18 18n ^5 / 

and the constant is given by 







01 9^01 r"+^ dr = — — / (pi r"+^ dr 
2 Jo 

For the associate eigenfunction in the geodesic ball B^{p) normalized to be positive 
and with L^-norm equal to 1, we have 



(14) e"/V(g) = 0i(y) 



where q is the point of M whose geodesic coordinates are e ?/ for y G -Bi , and G2 is defined 
implicitly as a solution of an ODE in [S]. Although we do not need its expression, for 
completeness we recall it: if we solve such ODE we found 

(15) G,M = _l.,.^,.M-c=^^^.(r). 



4. Known results 

Our aim is to perturbe the boundary of a small ball Bf{p) with a function v in order to 
obtained an extremal domain Bfj^^{p). The natural space for the function v is C^'"(S'"~^) 
but not all functions in this space are admissible because v must satisfy also the condition 

Volg5?+,(p)=Vol5i 

In order to have a space of admissible functions not depending on the point p, we use a 
result proved in [12] , that allows to use as space of admissible function the space 

C^"(5"-^) = S.V e C2'"(^"-i) : [ v = 

The result is the following: 

Proposition 4.1. (Pacard - Sichaldi [12]^' Let p G M . For all e small enough and all 
function v G C^°'{S^''~^) whose C'^'°' -norm is small enough there exist a unique positive 
function = (f){p,e,v) G C^'"(-Bf_,_^(p)), a constant A = A(p, e,t;) G M and a constant 
Vq = Vo{p, e, i;) G M such that 

Vo\gBf_,,{p)=Vo\Bi 
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where v := Vq + v and (f) is a solution to the problem 

Ag0 + A0 = in 5f+„(p) 
= on dBf^^ip) 



(16) 
normalized by 



[ 0Mvolg = l. 



In addition (j), A and vq depend smoothly on the function v and the parameter e and (f) = 
A = Ai and vq = when e = and v = 0. Moreover Vo{p, e, 0) = O(e^). 

Instead of working on a domain depending on the function v = Vq + v, it will be more 
convenient to work on a fixed domain Bi endowed with a metric depending on both e and 
the function v. This can be achieved by considering the parametrization of Bf^^{p) given 
by 



Y{y) 



where x is a cutoff function identically equal to when \y\ < 1/2 and identically equal to 
1 when \y\ > 3/4. Hence the coordinates we consider from now on are y G -Bi with the 
metric g := Y*g. 

Up to some multiplicative constant, the problem we want to solve can now be rewritten 
in the form 



;i7) 



Ag0 + A0 = in Bi 
</) = on dBi 



with 

ri8) / (^2^voi 



9 
Bi 



and 

(19) Volg(5i) = Vol5i 



When e = and = 0, a solution of f[T7j) is given by = </)i, A = Ai and vq = 0. In the 
general case, the relation between the function (p and the function is simply given by 

Y*(f) = and A = A . 
After canonical identification of dBl^,^,{p) with S*""^, we define the operator 
F(p, e, v) = g{V(j), u) Iq^s - ^ - [ giV^, u) dvol^^^ , 
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where v denotes the unit normal vector field to dBf_^_^, and (</>, vq) is the solution of ( IT6l) 
provided by the Proposition 14.11 Recall that v = vq + v. Schauder's estimates imply that 
F is well defined from a neighbourhood of M x {0} x {0} in M x [0, oo) x C^"(S'"~^) into 
C^"(S'"~^) (the space C^^"(S'"~^) is naturally the space of functions in (^^'"(S*""^) whose 
mean is 0). Our aim is to find {p,e,v) such that F{p,e,v) = 0. Observe that, with this 
condition, (p will be the solution to problem ([7]). 

We also have the alternative expression for F using the coordinates of Bi and the metric 

9- 

where this time is the the unit normal vector field to dBi using the metric g. 
For all V G C^°(S'"~^) let ip be the (unique) solution of 



(20) 



A^tp + Xii) = in Bi 
ip = — ci V on dBi 



which is L (i?i)-orthogonal to 0i, where Ci := 9^01 Define 

(21) H{v):= idr^ + C2v)\Q^^ 

where C2 = d'^(pi\r=i- We recall that the eigenvalues of the operator — A^n-i are given by 
/^i = J ('^ ~ 2 + j) for j G N, and we denote by Vj the eigenspace associated to fij. 

The following result shows that H is the linearization of F with respect to w at e = 
and V = 0: 

Proposition 4.2. (Pacard - Sicbaldi, [12] j The operator obtained by linearizing F with 
respect to v at e = and v = is 

H : Cf,;"(5"-i) ^ C,i^"(5""^) 

It is a self adjoint, first order elliptic operator. The kernel of H is given by Vi. Moreover 
there exists c > such that 

||w||c2,-(S"-i) < c||iJ(w)||ci."(5"-i) , 
provided w is L'^{S^~^)- orthogonal to Vq (BVi. 

Using the previous proposition, the implicit function theorem gives directy the following: 

Proposition 4.3. (Pacard - Sicbaldi, |12]j There exists eo > such that, for all e G [0, eo] 

and for all p G M, there exists a unique function v = v {p,e) G C^"(^"~^), orthogonal to 
Vo(BVi, and a vector a = a{p,e) G M" such that 

(22) F{p,e,v) + {a,-)=0 

The function v and the vector a depend smoothly on p and e and we have 



\a\ + ||w||c2,a(5n-i) < ce 



2 
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In other word, for every point p G M it is possible to perturbe the small ball Bflp) 
in a domain Bf_^_^{p), whose volume did not change, but with the (strong) property that 
F{p, e, v) (i.e. the Neumann data of its first eigenfunction minus its mean) is the restriction 
of a linear function (a, ■) on S"'~^. It is important to underline that this result does not 
depend on the geometry of the manifold, because it is true for every point p. 

Now, we have to find the good point p for which such linear function (a, ■) is the 
function. And in this research we will see the geometry of the manifold. 

5. Construction of small extremal domains 
For p G M, let us define the function 

'^eip) := A = \{p,e,v{p,e)) 
where A is given by f|T7|) taking v = v{p, e) given by Proposition 14.31 

Proposition 5.1. Fore small enough, the domain B^^^(^p^){p) is extremal if and only if p 
is a critical point of where v{p, e) = Vo{p, e, v{p, e)) + v{p, e). 

Proof. Recall that by definition 

F(p, e, v{p, e)) = giu, V0) - b 

where 
and 

/ ^dvolg.„ = 0. 

JdBi 

Moreover we know that 

F(p, e, v{p, e)) + (a(p, e), ■) = 0. 
In particular the domain B^^^f^^ ^^{p) is extremal if and only if a(p, e) = 0. 

Let us now compute the differential of \l/e. Let S G TpM and 

q := Expp(tS). 

For t small enough, the boundary of Bf_^_^^^ (q) can be written as a normal graph over the 
boundary of i?f_,_^(p (p) for some function /, depending on p,e,t and H, and smooth on t. 
This defines a vector field on dB^^_^^^ (p) by 



dt 



V 

t=0 



where u is the normal of dBf_^_^^^ ^^{p). Let X be the parallel transported of S from 
geodesic issued from p. As the metric g is close to the Euclidean one for e small, there 
exists a constant c such that for all e small enough and any H the estimation 
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= - / |9(W,«')l'9(2,!>)dvoli,.. 

JdBi 



holds. The variation of the first eigenvalue, see Proposition I2.H gives 

d 

We thus obtain 

(23) = - / [-{a{e,p),-) +hf g{Z,u)dYo\g,^ 

JdBi 

Recall that the variation we made is volume preserving, i.e. 



/ g{Z, u) dvolg;^ = . 

JdB^ 



Then it is easy to see that if a = then Dp^^ = 0. This proves one implication. 
For the reverse implication, assume now that Dp'^^ = 0. From f l2^ we have 

(24) 2b [ {a{e,p), ■) g{Z, z>) dvol^,^ = / (a(e,p), ■)^g{Z, z>) dvol^,^ 

JdBi JdBi 

for all H. It is easy to see that for all e small enough there exists a constant c such that 

(in fact the left hand side vanishes when e = 0, the metric g and the Euclidean one differ 
by terms of order and the normal vectors differ by terms of order e). Now we choose 
E = ba = b{p, e) a{p, e) and we get 

g{Z,u) = b{a,-) + eA 



where |y4| < c||6a||g. Using this equality in equation fl24|) . we deduce that for all e small 
enought there exists a constant C independent on e and a such that 



26^ 



/ (a, Y dvolg;,^ <C\b\{e \\af + \\af + e\\a\ 

JdBi 



Now the left hand side is bounded by below by 6^ so finally we obtain 

b^\\af <C\b\ {e\\a\\ + \\a\\ + e) \\a\\\ 
Observe that we cannot have 6 = because it would imply that 



A / dvolg = / dvolg = 6 = 

Bi Jb, 



but > on i?i so A = which is not possible because it is the first eigenvalue of 
the Laplace-Beltrami operator with Dirichlet boundary condition. As ||a|| = O(e^), 
then for e small (recall 6 7^ 0) we obtain that a = and this concludes the proof of the 
proposition. □ 
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We now define 

n [n + 2) + 2Ai 

Propositions 14.31 and 15.11 completes the proof of tlie first part of Tlieorem II. 1[ In tlie 
following sections, we will prove the second and the third parts of Theorem 11.11 and for 
this we have to find an expansion in power of e for ^'e(p). Such expansion will involve the 
geometry of the manifold. 



6. Expansion of the first eigenvalue on perturbations of small geodesic 

BALLS 

In this section we want to find an expansion of the first eigenvalue A = X{p, e, v) in power 
of e and v, where p is fixed in M. In a second time, we will use the function v = v{p, e) given 
by Proposition 14.31 in order to find an expansion of X{p,e,v{p,e)) in power of e. Keeping 
in mind that we will have v = 0{e^) we write formally 

A(p,e,t;) = A(p,0,0) + 9,A(p,0,0)e 

+diXip,0,0)v + ^df\ip,0,0) 
+d,di\{p, 0,0)ev + ]: df\{p, 0, 0) 

D 

+ ^ dfX{p, 0, 0) + ^ d%X{p, 0, 0) + ^ dt\{p. 0, 0) e' 

We thus study all of theses terms. 
Lemma 6.1. We have 

a,A(p,0,0) = 

\ dfXip, 0, 0) = -:|(l + 2 -j^) =: Ao 
2 6 \ n (n + 2) / 

df\{p, 0,0) = 

^dfXipAO) = A + Ai f^--— ^^"j =:A 
24 \ Un 36 [n + 2) J 

where the constants A and W are given in / flT]) and l[T^) . 

Proof. It suffices to find the expansion of X{p, e, 0) in power of e. First we have to expand 
Vq{p, e, 0) and this can be done by using expansion (fTOl) . keeping in mind the definition of 
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the metric g and the fact that when v = the constant vq is given by the relation 



We find 



where 



^0 

A 



VohBi = VolBi = — 
n 



Wo 

— ({n + 2)AoWo + ^—Alu^ + W 

UJn \ I 



Now we use expansion (fT^ replacing e by e(l + v^)- We obtain 

A = Ai + Aoe2 + Ae4 + O(e5) 

where 

Ao = Ao-2AiAo = -^ + 

6 \ n (n + 2) 



A 



A- Ai(2A-3A2) = A + Ai 



This concludes the proof of the result. 
Lemma 6.2. We have 

5,A(p,0,0) 
dfX{p,0,0){v,v) 



-2ci / vH{v) 



□ 



where H is the operator of Proposition whose expression is given by l[21\). and Ci : = 
<9r0ilr=i the constant defined in 

Proof. Let Qo = Bi be the unit ball of M", and let Vtt = -B(i+t,p+tt-), where we recall that 
fgn-i V = and Vq = Vo{t) is chosen in order that Vol fit = Volfio = We have 



and 



d,X{p,0,0)iv) 



dlX{p,0,0){v,v) 



dt 



dt^ 



Ar 



t=0 
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where Aqj is the first Dirichlet eigenvalue of ilf The expansion of Vol ilt directly prove 
that vq = (9(t^). In fact, in polar coordinates, we have 

Vol a = / / r""Mrd^ 



- / {l + vo + tv^de 

n Jgn-l 

- [ [(1 + VoT + nil + voT-^ t V + o{f)] de 



n 



n 

Differentiating this expression with respect to t, and keeping in mind that vo{0) = 0, we 
obtain that vo(t) = For y & Qq and t small, let 

h{t, y)=(^l + vo + t xiv) V (^^^ ^ y 

where x is a cutoff function identically equal to when \y\ < 1/2 and identically equal to 
1 when \y\ > 3/4, so that h{t,flQ) = Qf We will denote the t-derivative with a dot. Let 
V{t, h{t, y)) = h{t, y) be the first variation of the domain Vlt. Let v be the unit normal to 
dVlt and let a = {V, v) the normal variation about dVtt. Let A be the first eigenvalue and 
the first eigenfunction of the Dirichlet Laplacian over VLt normalized in order to have 
norm equal to 1. From Proposition 12.11 we have 



A = -/ {d^Ya 

JdQt 

where d^(f) = (V0, u). At t = and on the boundary, we have = 0i, d^(j) = dr4>i = Ci, 
a = V. Then A(0) = 0. This proves the first part of the Lemma. 

We can use now equality (l33l) of Proposition 110.11 of the Appendix (with / = (9^0)^ cr) 
in order to derivate this formula with respect to t. We obtain 



A 



/ \id,(PY {a + ad,(l) +H a^) + 2a (9,0 9,0 + a 8,^ dl 
Jdnt '- 



'dQt 

where H is the mean curvature of dQf. Now the second variation of the volume of Qf is 



Vo\nt= {& + ad^a + Ha') = 0. 

JdQt 

Such equation can be obtained differentiating equality (!32|) of Proposition 110. II with / = 1, 
using equality (l32l) of Proposition 110.11 with f = a. On the other hand, at t = and on 
the boundary, we have = 0i, 9,0 = 9^01 = Ci, 9^0 = 9^0i = C2 = — (n — l)ci, a = v, 
<p = tp, where ip solve fl20|) . We obtain 

A(0) = -2 ci / V {drip + C2v) = -2 ci / v H{y). 
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The proof of the Lemma follows at once. □ 
Lemma 6.3. We have 

d,diX{p,0,0) = 



6 J gn-l 



where G is the vector ofTpM whose geodesic coordinates are y E , according with 
and Ci := dr(t)i\j.^i is the constant defined in ^21 



In order to prove this lemma, we start with a preliminary result. The formulas for the 
geometric quantities we will consider are potentially complicated, and to keep notations 
short, we agree on the following: any expression of the form Lp{v) denotes a linear combi- 
nation of the function v together with its derivatives up to order 1, whose coefficients can 
depend on e and there exists a positive constant c independent on e G (0, 1) and on p such 
that 

\\Lp{v)\\c^,c(^s^-i) < c ||t;||c2,a(sn-i) ; 

similarly, given a G N, any expression of the form Q^p\v) denotes a nonlinear operator in 
the function v together with its derivatives up to order 1, whose coefficients can depend 
on e and there exists a positive constant c independent on e G (0, 1) and on p such that 

IIQp^H^l) - Qp"n^2)||ci.«(5"-i) < C (||t;i||c2-(5n-i) + ||t;2||c2.«(5"-i))" ^ \\V2 - t;i||c2."(5"-i) 

provided ||t;j||c72,c«(5n-i) < 1, for i = 1,2. 
Lemma 6.4. We have 

d,Vo{p, e, 0)(tl) = 7^ / Ric{Q, Q)v+ I [O(e^) + Lp{v)] 

where 6 is the vector ofTpM whose geodesic coordinates are y G S*"""*^, according with Ij^. 

Proof. The expansion in e and v for the volume of the perturbed geodesic ball B^^^_^^^{p) 
is given in the Appendix of [13] (the corresponding notations with respect to [13] are 
Be{i+v)(p) ~ Bp,e{,—v) and n = m + 1). We have: 



n 

(25) -/ ^'+^^ / / mcp{e,e)v 



f gn — 1 J gn — l U J S"^ 

+ / {0{e'^) + e'Lp{v) + e^Qf\v) + Qf{v)) 



5,1-1 

n-l 



where is the vector in TpM whose coordinates are y E S"' , and Wq, W are given by 
(fTT]) . Putting V = vq + V m. expansion fl2^ . where /^n-i v = Q and Vq is chosen in order 
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that the volume of Jp) is equal to the volume of B^, we obtain 



n — 1 



UJn I 1 



n — 1 



vq - - e 



Ricp{e,e) 



Sn-1 



2 1 



v e" 

n-l 6 



Ricp{e, e) V 

-1 



where 



R 

Ric = Ric g 

n 



is the traceless Ricci curvature. In order to compute the expansion of 

vo := d^voip, e, 0){v) = vo{p, e, sv). 

we derivate with respect to s, at s = 0, equality 

VoLM,^ , -Jp) = Vol 5, 

using the expansion above. Recall that we know vq^p, e, 0) = 0{e^). We find 

1 







/ mcp{e,e)v + [ (Oie^) + Lp{v)) 



Finally 

vo = ^e^ [ RiCp{Q. e)v+ [ {0{e') + Lp{v)) 
This completes the proof of the Lemma. 

We are now able to prove Lemma 16.31 
Proof. (Lemma \6. 3\) . We make a development up to power 2 in e, of the function 

d 



ds 



A(p, e, sv) . 



From Proposition 12. 1^ we have 
d 



s=0 



X{p, e, sv) = - I g{V, P){giy<P, i>))Mvol3,„ 



where the deformation in a neighborhood of dBi is given by 

^(■s, y) = {l + vo{e, sv) + sv)y 

and 



V{y) 



dh 
ds 



s=0 



dyVo{p,e,0){v) + V 



y- 
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In that formula, the term ^(V0, 0) is computed with s = or equivalently v = 0. From 
the definition of g and the expansion of the metric g, when -y = we have 

1 



9^j 



{l + vo{e,0))-' dr = {l + vo{e,0)) 



-1 y_ 
\y\ 



The expansion of e, 0) is almost known: it suffices to replace e by e(l + Vq) in formula 
([HD. We have 

= 01 + eV2 + O(e') 

where 



f2iy) 



R,.y'y^-^\y? 



n 



^ + R,G,i\y\). 



Using the notation i = g''°'g^^Rakbi we thus have on dBi 



k „.l 



d 



d 



dy' 



;i-i,o(e,0))-'ci+€2S,/2 + 0(€= 



where on the boundary 



drf2{y) 



Cl_ 

12 



R^.y'y'-^-^ 



Rp 
n 



+ RpG',il). 



Now we have to expand the measure on the boundary. This is classical and can be done 
directly from expansion (|8]). We have 



dvolgiJgs^ = {l+VoY 



1 - -ffiCp(e,e)e' + C>(e') 



dvoll Qn-l 



where is the vector in TpM whose coordinates are y G 5""~^ and dvol|5n-i is the Euclidean 
volume element induced on S"'~^. For the term dvVo{p,e,0){v) appearing in V we use 
Lemma [6.41 We have 

J2 



dfiVo{p,e, 0){v) 



Geo. 



n JS 



Ricp{e, e)v + [O(e^) + Lp{v)] 



We finally obtain 
d 
d^ 

where 



s=0 



X{p, e, sv) = Ce' RtCp{e, e)v + [O(e^) + Lp{v)] 



„2 2 

n — __1 _ 1 I _! 
~ 6 12 6 
The proof of the Lemma follows at once. 



□ 
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Summarizing the results of Lemmas 16. ![ 16.21 and 16.31 we obtain the following: 
Proposition 6.5. Let p G M , let e and v he small enough. Then: 

-ci / vH{v) - [ RiCp{Q,Q)v 

+ / [0{e') + L,{v) + (v) + (v)] 

where is the vector of TpM whose coordinates are x G 5*""^ according with and we 
agree with the convention about Lp{v), Qp (v) and Qp (v) we gave before. 

Proof. It suffices to put together the results of Lemmas 16.11 16.21 and 16.31 □ 

7. Localisation of the obtained extremal domains 
Now we want to find the expansion of the function \E'e(p) in power of e. Recall that 

^e(p) = A(p,e,w(p,e)) 

In order to find such expansion we will relate the first term in the expansion of v{p, e) to 
the curvature of the manifold at p. 

The first term of the expansion of v{p, e) is related to the traceless Ricci curvature at p, 
as stated by the following: 

Proposition 7.1. We have 

v{p, e) = --|^ Rzcpie, 0) e' + 0{e') = -—^ ^ ^zc,(0, 0) + 0{e') 

where is the vector of TpM whose geodesic coordinates are y G 5'"^"'^ , according with 
and a2 is the eigenvalue of the operator H defined in Proposition ^.^ associated to the 
eigenspace V2. 

Proof. Let us recall that 

F{p, e, v{p, e)) + {a{p, e), ■) = 0. 

where 

\\^{p,e)\\c2.^(^S"-^) + < ce^ 

Now, because F{p,e,0) = C(e^) and because v{p,e) = C(e^), we can write 

F{p,e,v) = F{p,0,0) + d,F{p,0,0)e 

+d,F{p, 0, 0) V + ]^dlF{p, 0, 0) + 0{e^) 
= H{v) + ^d',F{p,0,0)e' + O{e') 
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In the computation of the mixed derivatives of A in the proof of Lemma 16131 we have aheady 
computed the expansion of g(V(f), z>) for v = 0, so we directly deduce 

Rp 
n 



Then we have 

a,2F(p,O,O) = |ffiCp(0,0) 

Writing 

a = Ope^ + C(e^) 

and 

and considering the expansion of F, from equation fl22l) we obtain 

(26) H^vp) + ^RiCp{Q,Q) = -{ap,-) 

We know that v, and hence Vp, is L^-orthogonal to Vq © Vi (see Propositions 14.31) . Observe 
that i?ic(e, 6) is L2(^"-i)-orthogonal to Vi since the function 6 — )■ Ric{Q, 6) is invariant 
when is changed into —0 and hence its L^-projection over elements of the form g^E, 0) 
is for every S. Then Ric{Q, 0) is L2(S""~^)-orthogonal to Vq © Vi. In fact Ric{Q, 0) is 
the restriction on S*""^ of a homogeneous polynomial of degree 2 which has mean 0, and 
then it is an eigenfunction for — A^n-i with eigenvalue 2n. As H preserves the eigenspaces 
of — A571-1 and his kernel is given by Vi (see Proposition 14. 2p . we have that there exists a 
constant 02 7^ such that 

H (^Ric{e, 0)) = as R^c{e, 0) 

From fl26l) we obtain 

-{ap,.) = H (^vp + ^^Rzc{Q,e) 

i.e. (ttp, •) is in the image of H. But it belongs also to the kernel of H, and then Op = 
and 

(27) H(vp + -^Rzc{e,Q) 

Now we remark that i Vp -\ — Ric{Q, 0) ) is orthogonal to Vq © Vi, and then 

V 12a2 / 

(28) y^ = _^i?,c(0,0) 

12a2 

In order to complete the proof of the proposition we use equation (13T|) and Lemma 18.11 
of the Appendix. □ 
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Now we are able to give an expansion for the function "^eip) in power of e. 
Proposition 7.2. We have: 
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(29) 



itp 



where Aq is defined in Lemma WH (note that ^ is well defined also when Rp = 0), and the 
function r can be written as 

r = Ki WRiemf + K2 \\Ricf + K^R^ + AgR 

for some constants Ki only depending on n. 



Proof. Replacing v with its expansion given by Proposition 17. II in the expansion of A given 
by Proposition I6.5[ we obtain 



Ai + Ao + A - Ci 



V I H{v) + ^ fficp(e, 6)) + 0{e') 



5" 



Ai + Aoe^ + Ae^ 
Ai + Aoe'+Ae^ 
Ai+Aoe'+Ae^ + 



"1 ,4 



ei (ffic,(e,e))^ + o(6^ 

144^2 Js"-! 



1 



72 a2n{n + 2) 
Ai 



\Ric„ 



n 



-Rl]+ 0{e' 



\Ri Or, 



Rl + 0{e') 



36(n + 2)(n- Ai) v ' ^ 
where we used (l28ll from the second to the third line, the following two geometric formulas 

ffic(e,e) = —R„ 



n 



{Ric{Q,Q)f 



n(n + 2 



^ {2\\RiCp\\^ + R: 



■p) ' 



whose proofs can be found in [13], from the third to the fourth line, and the computation 
of a2 given in ( 13T]) and Lemma 18.11 to deduce the last line. Define 



Tp Rp Aq 



A + 



^ ^ ( llffi. f - -Rl 

36(n + 2)(n- Ai) V n " 



A + Ai 



2W 



Rl 



Ai 



Un 36^2 (r2 + 2)7 36(n + 2)(n- Ai) 
Recalling the definition of W and A given in (ITT]) and (fT3l) . we obtain that 
rp = i^i llffiempip + K2 \\RiCpf + K^RI + {AgR)p 



Ricpf --Rl 
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where 

K, = , . . _ ( 18c^+ 



n(n + 2) + 2Ai V 10(n + 4) 

1 /35 2 4Ai riAi 

= . ...... I — + 



(30) 



is:. 



n(n + 2) + 2Ai V 3 15(n + 4) 6(Ai - n) 

1 5n — 3 2 Ai Ai 



3 



n(n + 2) + 2Ai V 3n 6(n + 4) 6n 6(Ai - n) 



n(n + 2) + 2Ai V5 5(n + 4) 

and formula ( 129|) follows at once. The fact that the constants Ki depend only on n comes 
immediately from the computation of (? by Lemma 18.21 in the Appendix: 

2 _ (n + 2) [2Ai + r2(r2 - 4)] 

12 AiWn 

This completes the proof of the proposition. □ 

Remark 1. We remark that ii'i > in order to justify our discussion about critical point 
of ||/22em|| for Einstein metrics in the introduction. 

Now recalling that 

_i*,(p)-Ai 6n(n + 2) - Ai 



<l>(p,e) = i?pAo 



" ° e2 n(n + 2) + 2Ai 

the proof of the second and third part of Theorem 11.11 follows at once. 

8. Appendix I : On the first eigenfunction in the unit Euclidean ball 

In this Appendix we state and prove some relations between the first eigenfunction and 
the first eigenvalue of the Dirichlet Laplacian on the unit ball. 

Lemma 8.1. Let 

ci = </>'i(l) 

where x — )■ (/)i(|x|) is the first eigenfunction of the Dirichlet Laplacian on the unit ball, 
normalized in order to have L'^-norm equal to 1. Then 



/2Ai 

where Ai is the first eigenvalue of the Dirichlet Laplacian on the unit hall. 
Proof. Recall that (pi is the solution of 



n — 1 



■ 0; + Ai0i = 
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with normalization 

1= / ^l{\x\) dx = cOn [ (0i)2r"-Mr = / ^i^^r^dr 



Bi 



n 



and 

Ai= / |V0i(|a;|)|Mx = a;„ r(0;)V"-Mr 

J Bi Jo 

Now let us compute 

(r"(0;)2)' = nr"-i(0'i)2 + 2r"(/);0'/ 



= (2 - ra) r"-i - 2Ai r" (p[ (pi 
Integrating this relation between and 1 we obtain 

2 2Ai 
Ci = 

The proof of the Lemma follows at once, keeping in mind that Ci is negative. □ 
Lemma 8.2. Let 

2 Jo 

where x — )■ 0i(|x|) is the first eigenfunction of the Dirichlet Laplacian on the unit ball, 
normalized in order to have L'^-norm equal to 1. Then 

2 _ (n + 2) [2Ai + n(n - 4)] 

12 AiWn 

where Ai is the first eigenvalue of the Dirichelt Laplacian on the unit ball. 
Proof. We have 

= f 4>l r"+i dr = - f 01 0; r^+^ dr 

2 Jo Jo 

Recall also that 

71 — 1 

0;' + 0; + Ai0i = o 

r 

with 01 (1) = 0, and 0i is normalized by 

(ojnY^ = I 02r"-Mr = -- I 0i0'ir"dr 
Jo n Jo 
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We first compute 

= (n + 2) r"+i - 2r"+2 0; 0; + 



2r" 01 0; + r"+i (0;)2 - Ai r"+i 0^ 



= (4-n)r"+i(0;)2-2Air"+2 0;0i 
Integrating this relation between and 1 we find 

cl = {4-n) rr"+i(0;)2 + 2Aic' 
Jo 

where ci = 0'i(l). We now compute 

(r"+i 01 0;)' = (n + 1) r'^ 0i + r"+i (0;)^ + r"+i 0i 0'/ 

/r; — 1 

= (n + l) 01 0; + r"+i (0;)2 - r"+i 0i 0^ + Ai 0^ 



(r" 0?)' - n r""i 0^ + r"+i (0;)' - Ai r"+^ 0? 



Integrating this relation between and 1 we find 



Ai 



n + 2 



We thus have at the end 



n + 2 
12Ai 



c + 



n(n — 4) 



The proof of the Lemma follows at once from Lemma 18.11 



□ 



9. Appendix II: The second eigenvalue of the operator H 

Here we compute the eigenvalue 02 of the operator H associated to the eigenspace V2. 
When w is an homogeneous polynomial harmonic of degree 2 (abusively identified with its 
restriction to the unit sphere) we have Asn-iw = —^2W = —2nw and H{w) = a2W. We 
recall that 

H{w) = {drilj)\Qs^ +C2W= ((9^^)|gs^ -{n-l)ciw 
where ip is the solution of 

r AV' + Ai^ = in fii 

y ip = — ci w on dBi 
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which is L^(i?i)-orthogonal to (pi. Decomposing in spherical harmonics, we see that 
■ip{r,6) = h2{r) w{6) where 62 is the solution defined at of 

r'^b" + {n-l)rb' + {r"^ Xi-2n)b = in (0,1) 

6(1) = -ci = -<P[{1) 
From the definition of H, we see that 

a2 = b',{l) + 0;'(1) = 6^(1) + C2 = 6^(1) - (n - 1) ci 
so we have to compute 62(1)- Let us verify that 

&2(r) = -f^0i + -0; 
\ n r 



is the desired solution. Recall that 



thus 



Now 



and 



so 



Tl — 1 

^1 + — 0; + Ai0i = o, 

r 



-^1 Mil , ^ j,w \ , 2 ,„ 2 



^2 = - + 



,„ 77, — 1,, In — 1,, n— 11,, 1,„ 1,, 
^2 + & +A162 = —<P'i + ^0 +2-0'-2-0 

ly ly ly ^ ry ry ^ ry ^ ryrj 



ly ly ^ ly ry ^ ry ^ lyO 

1 fn-1 , , \ 1 , 
-2- +Ai0i -2-0 

fy-ij \ 'y j iy-J 

2n fXi 1 



r"^ \ n r 



2n , 
02 



And of course ^2(1) = ~Ci, so this is the desired solution. Finally we have 

62(1 = Ci 

n 



and 



^- Ai Ai -n /2Ai 

(31) q;2 = ci = W > 0. 

n n \ ojn 
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10. Appendix III : Differentiating with respect to the domain 

In this Appendix we recall a useful result that allows to derivate the integral of a function 
with respect to a parameter t that appears in the function and also in the domain of 
integration. The proof of the such a result can be found in [7], page 14. 

Proposition 10.1. Let Q a smooth bounded domain o/M" and 

h : (-r,r) xil^W 

a smooth function, where r is a positive constant, such that h{0,p) = p for all p & Q. Let 

/ : M X M" ^ M 

a smooth function. Let ilt = h{t,Q,o), V{t, h{t,p)) = ^{t,p) and N{t,q) the unit outward 
normal at q ^ dQf Then 

f = [ %dx+ f f{V,N)ds 
and 

(33) ^J^^fds = j^^ (^-^ + (V, N) (V./, N)+H (r, N) f^ ds 

where (■, ■) denote the scalar product in W^, s denote the area element of dQt o-nd H is the 
mean curvature ofdflf 
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